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A search for disappearance of active neutrinos over a baseline of 735 km was conducted using the NuMI neutrino 
beam and the MINOS detectors. The data analyzed correspond to an exposure of 3.18 X 10^'^ protons on target. 
The data are fitted to neutrino oscillation models in which mixing with one sterile neutrino is assumed. A 
comparison of the neutral-currcnt-likc spectrum at the FD with the expectation derived from the near detector 
measurement shows that the fraction of disappearing muon neutrinos converting to a sterile state is less than 
52% at the 90% confidence level. In addition, the possibility of decay of active neutrinos into sterile species 
occurring concurrently with neutrino oscillations was analyzed. Pure neutrino decay is disfavored at 5.4 ct as an 
alternate explanation to oscillations for the depletion of muon neutrinos at 735 km. The methodology employed 
in the analysis of neutral current neutrino events observed in MINOS is described and newly obtained results 
are presented. 



1. The NuMI Beam and the MINOS 
Detectors 

The MINOS (Main Injector Neutrino Oscillation 
Search) experiment is a complete long-baseline neu- 
trino oscillation study. The NuMI (Neutrinos from the 
Main Injector) neutrino beam created at Fermi Na- 
tional Accelerator Laboratory (Fermilab) is sampled 
first by the Near detector (ND), on-site at Fermilab, 
at 1 km from the target, and then by the Far detector 
(FD), 735 km away in the Soudan Underground Labo- 
ratory in Minnesota. Neutrino oscillation phenomena 
are studied by comparing the reconstructed neutrino 
energy spectra at the Near and Far locations. 

The NuMI neutrino beam is produced using 
120 GeV protons from the Main Injector. The protons 
are delivered in 10/is spills with up to 4.0x10^'^ pro- 
tons per spill. Positively charged particles produced 
by the proton beam in a graphite target (mainly 7r+ 
and K~^) are focused by two pulsed parabolic horns 
and are then allowed to decay in a 675 m long, 2 m 
diameter decay pipe. The target position relative to 
the first horn and the horn current are variable. The 
data employed in this analysis were obtained using 
the low energy beam configuration, in which the peak 
neutrino energy is 3.3 GeV and were recorded be- 
tween May 2005 and July 2007, corresponding to a 
total exposure of 3.18x10^° protons on target (POT). 
The charged current (CC) neutrino event yields at the 
ND are predicted to be 91.8% z/^ 
and 0.1% Ve. 

The MINOS detectors are designed to be as sim- 
ilar as possible in order to minimize systematic un- 
certainties. The detectors are fine-grained tracking 
calorimeters with an inch thick absorber layer of steel 
and a 1 cm active layer of plastic scintillator constitut- 
ing one "plane" . Each scintillator layer is constructed 
from 4.1 cm wide strips. Signals from the scintillator 
are collected via wavelength-shifting (WLS) fibers and 
carried by clear optical fibers to photomultiplier tubes 



6.9% V^, 1.2% ly. 



(PMTs). The ND and FD arc both magnetized with 
a current-carrying coil producing an average field of 
1.3 T in the fiducial volume. 

The 0.98 kton ND, 103 m underground, has 282 ir- 
regular 4x6 m^ octagonal planes. The detector con- 
sists of two sections, a calorimeter encompassing the 
upstream 121 planes and a spectrometer containing 
the downstream 161 planes. In both sections, one out 
of every five planes is fully covered with 96 scintillator 
strips attached to the steel plates. In the calorimeter 
section, the other four out of five planes are partially 
covered with 64 scintillator strips whereas in the spec- 
trometer section no scintillator is attached to the steel. 
The 5.4 kton FD, 705 m underground, has 484 octag- 
onal, 8 m wide instrumented planes. Due to the ND 
proximity to the target, the signal rate in the ND is 
^ 10^ times larger than in the Far detector. 



2. Neutrino Interactions in the MINOS 
Detectors 

There are two main types of neutrino interactions 
observed in the MINOS detectors: Charged Current 
(CC) interactions, proceeding through the exchange of 
a boson with creation of the associated charged 
lepton, typically defined by a long muon track ac- 
companied by a small hadronic shower at the event 
vertex caused by nuclear fragmentation; Neutral Cur- 
rent (NC) interactions, proceeding through exchange 
of a Z'^ with the neutrino leaving the detector, which 
appear as diffuse showers with typical length much 
shorter than the length of a muon track in a CC event. 
Correct identification of NC events in the MINOS de- 
tectors is difficult due to to short CC events with high 
hadronic inelasticity, featuring a short muon track of- 
ten concealed by the hadronic shower or easily con- 
fused with the short charged pion tracks also found in 
NC events. 
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Figure 1: Example of CC (top) and NC (bottom) inter- 
actions in the MINOS detectors. Longitudinal views UZ 
and VZ are shown in both cases. CC events are charac- 
terized by a long muon track and hadronic activity at the 
vertex, whereas NC events are shorter, displaying a diffuse 
hadronic shower. 



3. Near Detector Selection 

Events reconstructed in the fiducial volume of the 
Near and Far detectors are very similar. However, in 
the ND, an average of 16 neutrino interactions are pro- 
duced for each 1.8 /zs spill during typical running with 
intensities of 2.2 x 10^^ POT per beam spill. Recon- 
struction algorithms are designed to handle this high 
interaction rate, but for certain event subclasses short- 
falls have been identified and quantified using special 
studies including low intensity beam data and visual 
scanning. 

Reconstruction failures are classified into three 
main categories: i) split events, ii) leakage events, and 
iii) incomplete events. Split events occur when a sin- 
gle neutrino interaction results in two or more recon- 
structed events. Leakage events are due to incorrectly 
assigned event vertices causing neutrino interactions 
outside the fiducial volume to be reconstructed within 
it. The incomplete event category is a looser classifi- 
cation that refers to further types of failures in shower 
reconstruction to be described below. In all three cat- 
egories, the visible energy of a neutrino candidate may 
be underestimated, resulting in a background to NC 
events at low energies. As ND data are used to pre- 
dict the expected spectra at the FD, reconstruction 
failures specific to the ND must be minimized. A set 
of selection requirements based on spatial and timing 
criteria, as well as on activity on the ND sparsely in- 
strumented regions, was developed to reduce the 
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Figure 2: The distribution of the number of strips with 
non-zero pulse-height, per event, after all other selections 
are applied. A sizable contribution from poorly recon- 
structed events (hatched histogram) is observed at low 
strip counts. The event range accepted for the analysis 
is identified by the arrow. 



occurrence of these failure modes. The selection cri- 
teria can be summarized as follows: i) the time sep- 
aration between events, Ai, must exceed 40 ns; ii) if 
40ns < At < 120ns, the spatial separation between 
events, Az, must exceed 1 m; iii) the ratio between 
the number of active strips per event plane and total 
number of active planes in the event must be less than 
unity; iv) for events with less than 5 GeV of recon- 
structed energy in which the number of planes is larger 
in the reconstructed shower than in the reconstructed 
track, the number of event strips reconstructed in the 
detector's veto regions should be less than four, or 
else the total pulse-height in those regions must be 
less than 2MIP^; and v) the total number of strips 
reconstructed in the event must be more than four. 
Only events that satisfy these criteria are used for fur- 
ther analysis. After applying these requirements, the 
background of poorly reconstructed events having vis- 
ible energy below 1 GeV has been reduced from 34% 
to 8%. The distribution of the number of occupied 
strips per event is shown in Fig. [2l 



4. Neutral Current Event Classification 

After the selection criteria described above are ap- 
plied, the analysis proceeds by distinguishing NC 



^Minimum Ionizing Particle, equivalent to the detector re- 
sponse to a perpendicular 1 GeV-muon traversing one scintilla- 
tor plane. 
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events from CC events. The goal of the event clas- 
sification is to maximize the efhciency and purity of 
selected samples of NC and CC events, with efficiency 
defined as the number of true events of one type which 
are classified as that type, divided by the total number 
of true events of that type. Purity is defined as the 
ratio of the number of true events of one type selected 
to the total number of events selected as that type. 

A sample of candidate NC events is obtained by 
applying specific requirements on three classification 
variables: event length, expressed as the difference be- 
tween the first and last active plane in the event; num- 
ber of tracks reconstructed in the event; and track ex- 
tension, defined as the difference between track length 
and shower length. Events crossing fewer than 60 
planes and for which no track is reconstructed are 
classified as NC. Events crossing fewer than 60 planes 
that contain a track arc classified as NC if the track 
extends less than 5 planes beyond the shower. The 
values chosen maximize sensitivity for detection of 
sterile-neutrino admixture. Finally, events that are 
not classified as NC-like are checked to determine if 
they are CC-like according to classification procedures 
described in a previous MINOS publication These 
requirements are applied to both ND and FD to obtain 
NC and CC event samples. Distributions for the event 
length and track extension classification variables for 
data of the ND are shown in Figs. [3^ and[3jD. The data 
are plotted together with the prediction of the MINOS 
Monte Carlo simulation, which adequately reproduces 
the shapes of the classification- variable distributions. 
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Figure 3: Comparisons of ND data with Monte Carlo pre- 
dictions for distributions of the variables (a) event length 
and (b) track extension. Systematic uncertainties are dis- 
played as shaded bands on the Monte Carlo expectation. 
Events selected as NC-like are indicated by the arrows. 



5. Near to Far Extrapolation 

The predictions of the energy spectra of the NC 
and CC samples at the FD are based on the observed 
ND data and make use of the expected relationship 
between the neutrino fluxes at the two sites. The 
process of making the predictions is called "extrap- 
olation" and may be viewed as making corrections 
to the simulation of interactions in the FD based on 
the energy spectrum measured in the Near detector. 
The current analysis uses an extrapolation technique 
caUed the "Far over Near (F/N)" method This 
method makes the prediction of the FD spectrum by 
taking the product of two quantities: The ratio of the 
expected number of events from the Monte Carlo sim- 
ulation for each energy bin in the FD and ND spectra; 
and the number of observed ND data events. The F/N 
method prediction is robust against distortions arising 
from differences between data and Monte Carlo sim- 
ulation in the ND as these distortions are translated 
to the FD and do not affect the oscillation measure- 
ment [l|. For instance, for the case of the CC 
component of the NC and CC samples, the F/N ex- 
trapolation predicts the number of events at the FD 



for the i-th bin of reconstructed energy to be 



ppredict jiqdata 



(1) 



where _/Vf is the number of selected events in the 
i-th reconstructed energy bin in the ND and Nf^^ 
is the number of events expected in that bin from the 
ND Monte Carlo simulation. The F^*^*^ represents the 
number of events expected from the FD Monte Carlo 
simulation in the i-th bin of reconstructed energy and 
j-th bin of true neutrino energy. In the equation, Ej is 
the true neutrino energy and P^^^^^ the probability 
of muon neutrino transition to any other flavor. 



6. Systematic Uncertainties 

The main sources of systematic uncertainties in the 
analysis are: i) absolute scale of the hadronic energy; 
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ii) relative calibration of the hadronic energy in the two 
detectors; iii) relative normalization between the two 
detectors; iv) CC background in selected NC events; 
and v) uncertainties due to the ND selection require- 
ments in the ND event counts. Monte Carlo studies 
have been performed where for each single uncertainty 
the Monte Carlo spectrum is varied by ±1 standard 
deviation independently, in order to estimate the ef- 
fect of each on the extrapolated spectrum. Beam and 
cross-section uncertainties that are common to the two 
detectors effectively cancel when using the F/N ex- 
trapolation. 

The absolute hadronic energy scale has an uncer- 
tainty of 12%. This value is a combination of the un- 
certainty in the hadronization model and intranuclear 
effects (10%) and uncertainty of the detector response 
to single hadrons (6%). The relative energy scale be- 
tween the two detectors has an uncertainty of 3%. 
The relative normalization between the two detectors 
has an uncertainty of 4%. This is a combination of 
the uncertainties due to fiducial mass, live time, and 
reconstruction differences between the two detectors. 

To evaluate the uncertainties due to the ND selec- 
tion, the requirement that the total number of recon- 
structed strips in an event is at least four was shifted 
by ±1 strip. The effects on the reconstructed en- 
ergy spectrum were determined for each shift. The 
uncertainty has been estimated to be 15.2% for 
£;reco <0.5GcV; 2.9% for 0.5 GeV < i^rcco < 1.0 GcV; 
0.4% for 1.0GeV< E.^^o < 1.5 GeV and is negligible 
fpr higher visible energies. 

Finally, the uncertainty in the number of CC back- 
ground events is determined using ND data taken in 
several different beam configurations. Using the ob- 
served differences in energy spectrum between the low- 
energy beam configuration and each of the other beam 
configurations along with information from Monte 
Carlo simulation of each configuration, the uncer- 
tainty on the CC background number is found to be 
15% for all energy ranges at the Near and Far detec- 
tors 

7. Search for Sterile Neutrinos 

Results from the MiniBooNE experiment Q , which 
looked for z^^ — > z^g or — > Ve appearance on a 
short baseline of 540 m using a neutrino beam with 
a mean energy of ^ 700 MeV, mostly rule out the 
potential scenario of oscillations at the mass-squared 
1 eV^ scale, a possibility reported by the LSND ex- 
periment @ . Scenarios including active-sterile neu- 
trino mixing with one or more sterile neutrinos have 
been put forward to explain the LSND result. Such 
scenarios, for which sterile neutrino oscillations occur 
over short baselines, are now essentially ruled out as 
an explanation for LSND 0], but the possibility of 
sterile mixing over long baselines remains. MINOS 



searches for sterile neutrino oscillations in the atmo- 
spheric sector, in a long baseline setting, by studying 
disappearance of NC events measured at the FD rel- 
ative to the flux observed at the ND. NC events are 
not affected by — > Vr oscillations, but would suffer 
an energy-dependent depletion if i/^ Vg oscillations 
were to occur. 

7.1. Active Neutrino Disappearance 

A three neutrino analysis assuming oscillations oc- 
cur only among active flavors is carried out to search 
for active neutrino disappearance. A total of 388 data 
events are selected as NC in the Far detector. The 
measured and predicted i?rcco spectra at the FD are 
shown in Fig. |4l The prediction assumes oscillations 
with the values of IAm32| and 623 previously mea- 
sured by MINOS Q. Although this analysis is not 
capable of isolating an electron neutrino appearance 
signal, it must take i/^ Ve oscillations into account 
because the classification criteria of this analysis in- 
clude Ve CC interactions in the NC enriched sample 
with nearly 100% efficiency. This is done by compar- 
ing the observed NC spectrum to two predictions, one 
that assumes null Vf, appearance, and another that 
assumes an upper limit for the Ve appearance rate in 
the FD calculated with the normal neutrino-mass hi- 
erarchy, 6113 = 12°, and 5 ~ 37r/2. The choice of 6'i3 
corresponds to the 90% confidence level upper limit 
established by the CHOOZ reactor experiment for 
the |Am32| value measured by MINOS Q. As seen in 
Fig. m the observed spectrum matches the prediction 
based on oscillations among the three active flavors 
quite well over the full range of allowed values of ^13 . 

The agreement between the observed and predicted 
NC spectra is quantified using a statistic, R: 

^= ? > (2j 

where, within a given energy range, A'oata is the ob- 
served event count, Bqq is the extrapolated CC back- 
ground from all fiavors, and S'nc is the extrapolated 
number of NC interactions [3] . Because the disappear- 
ance of occurs mainly for true neutrino energies 
< 6 GeV the data are separated into two sam- 
ples. Events with i?reco < 3 GeV are grouped into a 
low-energy sample while events with 3 GeV < -Ereco < 
120 GeV arc grouped into a high-energy sample. The 
values of S'nc a-nd contributions to -Bcc for the two en- 
ergy ranges are shown in TableU The measured values 
of R for each energy range indicate that neutrino oscil- 
lations among the active flavors describe well the ob- 
served data. Over the full energy range, — 120 GeV, 
a value of i? = 1.04±0.08(stat.)±0.07(syst.)-0.10(z^e) 
is measured, corresponding to a depletion of the to- 
tal NC event rate assuming null (maximally-allowed) 
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Figure 4: The reconstructed energy spectrum of NC se- 
lected events at the FD (points with statistical uncer- 
tainties). The Monte Carlo prediction assuming standard 
three-flavor oscillations is also shown, both with (dashed 
line) and without (solid line) Ve appearance at the CHOOZ 
limit. The shaded region indicates the 1 standard de- 
viation systematic uncertainty on the prediction. The 
hatched region shows the Monte Carlo prediction for the 
background of misidentified CC events in this sample. 



Table 1 Values of R, A^Data, S'nc, and the contributions to 
Bcc for various reconstructed energy ranges. The num- 
bers in parentheses are calculated including Ve appearance 
at the CHOOZ limit, as discussed in the text. The first 
uncertainty in the value of R shown is the statistical uncer- 
tainty, the second is systematic uncertainty, and the third 
is due to possible appearance. 



-Ercco (GeV) A'Data SnC -^CC 




^cc 


0- 


3 


141 


125.1 13.3 


1.4 


2.3 (12.4) 


3- 


120 


247 


130.4 84.0 


4.9 


16.0 (32.8) 


0- 


3 


R = 


0.99 ±0.09 ±0.07 


-0.08 




3- 


120 


R = 


1.09 ±0.12 ±0.10 


-0.13 




0- 


120 


R = 


1.04 ±0.08 ±0.07 


-0.10 





Ve appearance of less than 8% (18%) at 90% confi- 
dence level. 



7.2. Active-Sterile Neutrino Mixing 

The description of mixing between three active neu- 
trino flavors and one sterile neutrino requires the ad- 
dition of one mass eigcnstate, The choice of parame- 
terization for the expanded mixing matrix follows con- 
siderations in Ref. Two cases are considered for 
the placement of 1^4 in the neutrino mass spectrum: i) 
TO4 = mi; and ii) 7714 ^ 7713. In the first case, the first 
and fourth mass eigenstate are assumed to be degen- 
erate, as are the second and fourth mass eigenstates. 
These degeneracies imply that 6'i4 = 6*24 = 0° and the 



corresponding oscillation probabilities are given by: 



(3) 



where a = e,r, or s and A31 = (m| — m\)L/{AE). 
In the m4 3> case, the fourth mass eigenstate is 
assumed to be much larger than the third. Conse- 
quently the values of sin^ A41 and sin^ A43 average to 
i. Additionally, sin2A4i and sin2A43 average to 0. 
In this model Am43 is assumed to be O(eV^) such that 
the regime of rapid oscillations and thus the averages 
mentioned are valid at the Far site, while ensuring no 
observable depletion of occurs at the Near detec- 
tor. Using the above simplifications, the oscillation 
probabilities are: 



P. 



1-4| 1/7^3 



i7^4r sin^ A31 



(l-IC/^r) : 



= 47^ 



\UMU^^\' + U;,U^,U^3U:3 ) sin^ A31 



The data are compared to Monte Carlo predictions 
based on the probabilities in Eqs. ([3|) and ([4]) using a 
statistic including nuisance parameters for the five 
systematic uncertainties described above. The best fit 
values for the mixing angles in the two models as well 
as the for each are shown in Table HIl 



Table II Best fit points and uncertainty ranges obtained 
for the active- sterile oscillation models. Results are shown 
with and without appearance at the CHOOZ limit. All 
angles are given in degrees. 



(4) 



Model 


613 


xVd.o.f. 


O23 


O24 


ds4 


7714 — nil 



12 


47.5/39 
46.2/39 


45.0t«;° 
47.li?i% 




23.0t^l^ 


7714 > "73 



12 


47.5/38 
46.2/38 


45.0l«3:° 


u-"-o.o 

o+''-^ 

'-'■"-0.0 


^+^^■'' 
'J-J^-o.i 

23.0124.? 



From these results, a limit can be set on the mix- 
ing angles, ^34 < 38° (56°) at 90% confidence level 
for the 777,4 ~ fTT-i model. The number in parentheses 
represents the 90% C.L. limit obtained when maximal 
Ve appearance is allowed. For the 777.4 3> 7713 model, 
6I24 < 10° (11°) and 6*34 < 38° (56°) at the 90% confi- 
dence level. 

A perhaps more straightforward way to quantify the 
coupling between the active and sterile neutrinos is to 
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Table III Best fit points and uncertainty ranges obtained 
for tiie relevant parameters of the oscillation with decay 
model. The result obtained for the pure decay scenario, 
Am^2 ~* 0, is also presented. Angles are shown in degrees. 



Model 


xVd.o.f. 


a (GeV/km) 


e 


Osc. with Decay 


47.5/39 


0.001^0° X 10"' 


45.0i^»,f 


Pure Decay 


76.4/40 


4.6tl l X 10"' 





respect to previous MINOS analyses of neutrino de- 
cay, which used the CC FD spectrum only, since any 
neutrino decay would also deplete the NC spectrum. 
A Ax^ of 28.9 is obtained for the pure decay scenario. 
Thus, a pure neutrino decay model with null oscilla- 
tions, as considered in Ref. [8;], is disfavored at the 
level of 5.4 standard deviations, an improvement of 
1.7 standard deviations on the previously published 
value. 



determine the fraction of disappearing j/^ that transi- 
tion to i^s- That fraction is expressed as 



fs 



1 - p. 



(5) 



The 90% confidence level limit for fs is determined by 
selecting a large number of test values of 6*23 and 6*34 
from Gaussian distributions with mean and a given in 
Table El The value of fs that is larger than 90% of 
the test cases represents the limit. For both models, 
the value corresponding to the 90% confidence level 
is fs < 0.52(0.55), with the value in parentheses in- 
dicating the value obtained for maximally-allowed Ve 
appearance in the beam. Therefore, approximately 
50% of the disappearing can convert to Vg at 90% 
confidence level as long as the amount of i>i^ appear- 
ance is less than the limit presented by the CHOOZ 
collaboration. 



9. Outlook 

The analysis described here utilizes a exposure of 
3.18x10^° POT. However, MINOS has aheady accu- 
mulated more than 7x 10^° POT of NuMI beam data. 
An analysis of neutral current events in the additional 
data is underway with expected improvements of the 
limits on mixing angles and of the sensitivity to sterile 
fraction. Updated results are expected to be available 
in early 2010. 
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8. Neutrino Oscillations with Decay 

Neutrino decay, as an alternative or companion pro- 
cess to neutrino oscillations offers some capability for 
reproducing neutrino di sapp earance trends [ll|. The 
model investigated here |12l | includes neutrino oscilla- 
tions occurring in parallel with neutrino decay. In this 
model, the survival and decay probabilities are: 



P 



tj.li 



cos 



sm 
2 



decay 



2 cos W sm 

1 - e 



'"3 J- 

e ^^3^ cos 



sin^ d. 



'32-' 



2E 



(6) 



where T3 is the lifetime of the 1^3 mass state and 9 is the 
mixing angle governing oscillations between and Vt- 
The limits r3 ^ 00 and Ato^j ~^ correspond to a 
pure oscillations or a pure decay scenario, respectively. 
The best fit values extracted for 6 and the parameter 
a = m^/r^ using this model are summarized in Ta- 
ble [Till The results are consistent with maximal mix- 
ing {9 = 45°) and with no neutrino decay {a = 0). 
The 90% C.L. limit found for the neutrino decay life- 
time is r3/TO3 > 2.1 X 10^^^ s/cV. 

Both the NC and CC FD spectra are included in 
the fits, therefore additional sensitivity is gained with 
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